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Using limits on photon flux from Dwarf Spheroidal galaxies, we place bounds on the parameter 
space of models in which Dark Matter annihilates into multiple final state particle pair channels. 
We derive constraints on effective operator models with Dark Matter couplings to third generation 
fermions and to pairs of Standard Model vector bosons. We present limits in various slices of 
model parameter space along with estimations of the region of maximal validity of the effective 
operator approach for indirect detection. We visualize our bounds for models with multiple final 
state annihilations by projecting parameter space constraints onto triangles, a technique familiar 
from collider physics; and we compare our bounds to collider limits on equivalent models. 


INTRODUCTION 

In this era constraints on Dark Matter models are being synthesized from multiple experiments. There has been 
much recent work in collider physics, focusing on Dark Matter (DM) models, which include both UV complete theories 
and Effective Field Theory (EFT) scenarios. The same models studied in collider physics imply detectable signatures 
from Dark Matter annihilation in space. Due to gauge invariance or other theoretical considerations, many of these 
models, both EFTs and simplified models, predict couplings between Dark Matter and multiple species of Standard 
Model particles. Thus, Dark Matter may be produced in many correlated final state channels at colliders, and may 
have multiple final state annihilation channels in space, which would contribute to total detectable photon or positron 
flux for satellite experiments. 

In this work we explore the indirect detection bounds from Fermi-LAT dwarf spheroidal galaxies |I] on models 
where Dark Matter annihilates into multiple final state channels. These bounds are among the tightest constraints 
on DM models. We study several EFT models with dimensions 6 and 7 effective operators. We choose operators 
which lead to unsuppressed DM annihilation rates in indirect detection processes, and which are being simultaneously 
studied in DM production processes at LHC. The dimension 6 operators we study are those which couple Dark Matter 
to third generation fermions pairs. The dimension 7 operators we consider are vector boson portals where there DM 
couples to multiple pairs of SM gauge bosons. 

The use of effective operators allows a great degree of model independence for Dark Matter studies, while capturing 
some of the important kinematic features of Dark Matter processes [HU. For some models which are completed by 
loops, EFT based calculations have so far provided the best means for study at colliders. The limits of the effective 
operator paradigm are becoming more clear for collider analyses. In particular, UV completions of models with low 
scale messenger portals are less probe-able by colliders, and models with low scale effective operator cut-offs may 
not be sensible at collider energies IMS]. However, we expect that EFT analyses are reliable at the scale of indirect 
detection where the center of mass energy of the annihilation process is the same order as the dark matter particle 
mass itself. We expect that the EFT treatment is valid down to much lower scales, perhaps for mediator sectors in 
range of 10’s of GeV, as opposed to colliders where mediators of some hundred GeV to just under a TeV may not be 
appropriate. 

We will visualize the bounds we set in two slices of the total parameter space, the plane of fixed DM annihilation 
rate, and the plane of effect operator coefficients where the total annihilation rate varies. For the regions of fixed 
annihilation rates, we will use the constraints to produce a 2-D visualization of the bounds on a triangle, a technique 
familiar from collider physics [9]. We also compare bounds set with the dwarf limits to those set by collider constraints, 
and discuss the validity limits of EFTs for both cases. 

The format of this paper is as follows, in Section I we will discuss dwarf constraints on photon flux from dark 
matter annihilations. In Section II we will analyze constraints on models with non-interfering final state annihilations 
and present triangular visualizations of parameter space. In Section III we will analyze a popular set of vector boson 
portal models with interfering final state channels. Section IV presents results along with collider constraints and 
discusses EFT validity. Section V concludes. 


INDIRECT DETECTION FROM DWARF SPHEROIDAL GALAXIES 

Dwarf spheroidal galaxies provide some of the tightest constraints on photon flux originating from dark matter 
annihilation as they are believed to contain a substantial dark matter component [TOl [IT] . This combined with their 
low astrophysical background makes them a good laboratory to search for dark matter. As no significant excess in the 
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photon spectrum has been observed from dwarf data, we use upper bounds on photon flux obtained from Fermi-LAT 
data pQ to place constraints on dark matter mass and couplings in the scenarios discussed in the introduction. 

The photon flux (photons cm -2 s -1 ) at the earth expected from annihilation of dark matter in the area of interest 
is given by 
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Where the J-factor (GeV 2 cm 5 ) is the line of sight integral of the dark matter density p, integrated over a solid angle, 

An 

J = f f p 2 (r)dldn'. (2) 
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The additional terms in Eq. [l] are dependent on the particle properties of the dark matter. Here, dNf/dE is the 
differential photon energy spectrum per annihilation to final state f, (av) / is the thermally averaged DM annihilation 
cross section and m x is the dark matter mass. 

The Fermi-LAT collaboration has presented dark matter constraints from observation of 25 dwarf spheroidal galaxies 
PQ. Of these 25 Milky Way satellites, 15 were used to place constraints on annihilation cross sections of dark matter 
particles with masses between 2GeV and lOTeV into various standard model channels. These annihilation constraints 
were derived assuming that the dark matter annihilates primarily to a single standard model channel. This is not 
realistic when considering most UV completions of Dark Matter models, even ones that are quite simple. Here we 
discuss obtaining bounds on dark matter mass and annihilation cross sections assuming annihilation into multiple 
channels. 

The bin by bin integrated y-ray energy flux upper limits at 95% CL for each of the 25 dwarf galaxies was presented 
in pQ. The J-factor for dwarfs is determined through dynamical modeling of their star densities and velocity dispersion 
profiles; we use J-factors presented in pQ. For a conservative limit we perform a very simple combination of six dwarfs 
with the highest J-factors. The limit in each bin is obtained from the dwarf which gives the strongest 95% CL upper 
limit on y-ray flux scaled by the J-factor. The combined limits are shown in Fig. [l] 
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FIG. 1: Combined 95% upper bounds on photon energy flux scaled by J-factor from 6 of largest J-factor dwarf spheroidal 
galaxies. 


We note that the Fermi-LAT collaboration has obtained its own combined dwarf limits on photon flux from DM 
annihilation. However, these limits assumed DM annihilation into one and only one visible channel at a time. Since 
our goal in this work is to analyze models that may have multiple final state annihilations we have forgone use of 
Fermi’s combination, though in future analyses, stronger bounds than ours may be obtained in each bin by ’’stacking” 
the total signal from all dwarves and subtracting a summed astrophysical background. 

We calculate the expected (av) / for each kinematically accessible annihilation channel, and obtain the differential 
gamma ray spectra for annihilation channels from the Mathematica code PPPC 4 DM ID [T2l H3] . As an example, we 
show the spectrum for several relevant SM pair annihilation channels for a lOOGeV mass dark matter particle in Fig. 
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FIG. 2: Gamma ray spectrum from DM annihilation into various channels for a lOOGeV DM fermion. 


[ 2 ] below. Here, the dimensionless parameter x is the gamma ray energy scaled by the dark matter mass. Analysis of 
many models requires calculation of differential gamma ray spectra from particles produced in non-symmetric pairs, 
for example, in the analysis of models of section IV we have calculated the differential spectrum of the Z 7 final state. 
Limits are obtained by binned comparison of total expected flux from all DM annihilations to our combined upper 
bounds. 


MODELS WITH INDEPENDENT ANNIHILATION CHANNELS 

We will now discuss the method of constraining the parameter space of models with multiple independent annihi¬ 
lation channels. We will first assume an effective Lagrangian which is the sum of several independent operators, each 
of which couples Dark Matter to one and only one pair of SM particles. Therefore the Lagrangian has the form 

C i = H i O i = H~xxX i X i . (3) 

Here the x is the Dark Matter, and X is some Standard Model particle with particle index i. Any specific operator 
will be gauge and Lorentz invariant and will have coefficient k/ A™, where the effective cut-off A appears with the 
appropriate power to make the operator dimension 4. 

The full parameter space of the EFT consists of the Dark Matter mass m x and the i operator coefficients ft/A™. 
Each point in this parameter space specifies a total DM annihilation rate, the specific ratios of the i final state 
annihilation channels, and the resultant 7 -ray flux. There are various slices of parameter space which can be studied. 
The first one we will consider is slices of the parameter space on which the total DM annihilation rate is held constant. 
Below we will show that along planes of parameter space with fixed annihilation rate, we will rule out masses and 
effective cut-offs below certain scales. 


Fixed Annihilation Rate 

The total annihilation rate, (crv ) tot , is simply a linear sum of the thermally averaged annihilation cross sections 
(■ crv)o i to particle X{ due to operator Oi, 

{ov) tot = N(av) Th = (av) 0l + {crv)o 2 + •■• (4) 

We will first fix the desired total annihilation rate. This rate may be anything we like, for simplicity we will 
consider it some numerical factor times the thermal annihilation rate N(crv}Th- This constraint drops us 1 dimension 
in parameter space; certain coefficient values k/ A™ will satisfy the constraint for any specific DM mass. Having 
fixed the total annihilations rate, we may then determine the limits of the operator coefficients which saturate the 
Fermi-LAT photon-flux bounds for any given Dark Matter mass. 

A natural choice for the total annihilation rate is the thermal rate. However, in the spirit of model independence, 
we will show Fermi-bounds on models with various annihilation rates which will correspond to models with various 
non-thermal histories. The complete theory will have to account for this history, as well as the presence (and absence) 
of specific operator coefficients. We note that models where the total visible annihilation rate is below the thermal 
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rate can easily be saved from the prospect of over-closure by the addition of invisible DM annihilation channels which 
restore the total annihilation rate to thermal. We discuss such models in the next subsection. 

The rate constraint confines us to a hyper-surface in the full parameter space, however we may greatly simplify the 
form of this constraint by dividing out by the total rate (crv} tot . We define the fractional annihilation rate as in 
to be Ri = (crv)i/ (av) t 0 t, giving a constraint which is linear in the partial rates Ri. 


R\ + i?2 T R 3 + • • • — 1- (5) 

The natural visualization for this parameter space of models with 3,4,5 annihilation channels etc. are a triangle, 
tetrahedron, 4-simplex, etc. 

As an example and for clarity of discussion we will analyze a simple model where DM may annihilate to multiple 
independent channels. We choose an EFT model where dark matter couples to the third generation of fermions via 
dimension 6 operators; 


£f = §* r * tn + Xf :xTxbTb + ^ 2 xTxrTr + ^X^v. (6) 

Here T specifies the whether the fermionic currents are scalar, pseudo-scalar, vector etc. Under the conventions 
in [3] they correspond to operators Dl, D2, etc. Annihilation rates will depend greatly on the Lorentz structure of 
the operator with operators leading to velocity and/or helicity suppression of the rate [15]. The strongest constraints 
from photon-flux measurements will apply to models with operators corresponding to unsuppressed DM annihilation 
rates. We have here included an annihilation into an invisible channel, specifically neutrinos. In setting constraints for 
fractional annihilation rates, we may consider this channel a stand-in for annihilation into any kinematically accessible 
invisible channel. The option of an invisible channel allows a total visible annihilation rate below the thermal rate 
while still avoiding over-closure of the universe. Recent work, for example os, has shown that these annihilations 
may be significant in some scenarios. The relevant annihilation diagram are shown in Fig. [3] 
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FIG. 3: Dark matter annihilation into pairs of fermions. 

Let us first assume that the dark matter is light and that annihilation into top pairs in kinematically forbidden. 
Here, the total annihilation rate is the sum of the thermally averaged annihilation cross section into b’s, r’s, and the 
invisible channel, 


(<™)tot = (<rv)b + {vv) r + (av) v . (7) 

The crucial point here is that each operator in Eq. [6] corresponds to only one channel, thus the constraint above has 
the form oc a (/^/A^) + b (ft r /A^) + c{k v ) A^) (where a, b, and c depend on dark matter mass and final state 
particle mass) with no terms of the form for i ^ j. Therefore, a triangle seems the natural choice for visualization 
of constraints where the sides of the triangle correspond to the Rb,R T ,R v axis. 

We see from Eq. [l] that once the annihilation rates are fixed, upper bound on y-ray flux corresponds to a 

lower bound on dark matter mass. To determine this mass bound we set the total annihilation cross section to a 
multiple of the thermal relic, (cry)to t = A"(cr?;)thm 5 then for a given point in (i2&, R r , R u ) space find the lowest mass 
which gives the maximum allowed y-ray flux determined by Fig. [l] The mass limits depend on the total annihilation 
rate which we have fixed; as we increase the total annihilation rate, the mass bounds become stronger. 

The results are given in the Fig. [4] for two total annihilation rates, the thermal rate and 10 times thermal rate. 
The blank spaces in the triangle correspond to a lower bound of 0 GeV or no lower bound. We see that the strongest 
bounds come from annihilation into b’s while the weakest bounds come from annihilation into invisible particles as 
expected. 
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FIG. 4: Lower bound on dark matter mass in GeV for annihilations into b quarks, r’s, and invisible particles for (av ) to t = {crv )Th 
on left and {crv )tot = 10(cru)Th on right. 

Here the points that form the triangle each correspond to a unique ratio of the partial annihilation rates Ri. At 
each vertex of the triangle, a single Ri = 1, that is, one channel saturates that total annihilation rate. Along the 
edges of the triangle only two annihilation channels contribute to the total annihilation rate, with the third partial 
rate set to 0. The form of Eqn. [6] was satisfactory to derive a lower DM mass bound for a fixed total annihilation 
rate. We will now pick a specific Lorentz structure for our operators and calculate bounds on the effective operator 
coefficients. As an example, we will choose to constrain the vector current operators D5, whose annihilation rate is 
neither velocity nor chirality suppressed. 

£f = %xYxh,J (8) 

A / 

Each mass that satisfies the rate constraint will correspond to specific values of the effective operator coefficient. 
With total annihilation rate fixed, points of lower mass will require lower effective cut-offs to satisfy the bound. We 
also plot on the triangle the lowest effective cut-off (largest operator coefficient) that satisfies the bound of total flux. 
Below we show the lower limits operator coefficients A* = Ai which saturate the bounds on photon flux. 



FIG. 5: Lower bound on operator coefficients A£ and A* for {av ) to t = {crv )Th on left and {av ) to t — 10(cru)Th on right. 

For the thermal total annihilation rate, we see that DM annihilating into mostly b’s or r's, must be above 5-15 
GeV in mass while the effective operator cut-off scales may not be smaller than a few hundred GeV. We point out 
that in our exclusion region, the effective operator paradigm appears to be justified as the effective cut-offs are well 
above the center of mass energy for the annihilation process. We show similar triangles for a total annihilation rate 
of 10 times the thermal rate in the appendix. 
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Constraints for 4 Independent Channels 

Next we consider the case of allowing more than three annihilation channels. We now allow the dark matter mass 
to be larger than the top mass opening up an additional annihilation channel. One can visualize the mass constraints 
on a tetrahedron where the “vertical” axis is R t , see Fig. [6j Each slice of the tetrahedron would correspond to a flat 
triangle with constraint equations given by 

(crv)tot - (<n;)t = (<™)b + (<™) r + (crv) v 

(1 — Rt) = Rb + R T + Rw (9) 



FIG. 6: Axis labeling for tetrahedron presentation of DM mass constraints when a fourth annihilation channel is allowed. The 
base corresponds to the top left of Fig. [7] 

In particular, the base of the tetrahedron would correspond to R t = 0 and the constraints shown in the top left of 
Fig. 0 Below we show lower bounds on dark matter mass for triangle slices where the fractional annihilation rate 
to tops is R t = 0, .3 and, .7. The blank areas in the triangles of Fig. [7] correspond to the region of parameter space 
where the kinematic bound is stronger then that obtained from indirect searches. 
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FIG. 7: The contours correspond to lower bounds (in GeV) on DM mass from annihilations into pairs of fermions (tt, bb , rr, isis). 
The total annihilation rate is set to 33 times the thermal relic for all three plots. The left image corresponds to Rt — 0, the 
right to R t — .3 and the bottom to R t = .7. 

For models with each operator corresponding to independent annihilation channels, this procedure may be iterated 
for any number of final state annihilation channels. In the next section we will discuss bounds on the parameter 
space of models where there is not a one to one correspondence between operators and annihilation channels. We first 
discuss the bounds on parameter space visualized on surfaces of fixed total annihilation rate. We will then discuss 
the bounds on surfaces of fixed operator coefficients. 


CONSTRAINTS IN MODELS WITH INTERFERING CHANNELS 

Many examples of models exist in which a single parameter or operator coefficient in the Lagrangian demands Dark 
Matter annihilation to multiple final state channels. Further there may be multiple operators or parameters in the 
Lagrangian which contribute to annihilations into the same channel. We will now analyze a model in which both 
of the previous statements are true. We consider a popular set of effective operator models for collider and indirect 
























studies (see, for example HMD), the vector boson portal models . We will consider a set of dimension 7 operators 
where a pseudo-scalar Dark Matter current couples to the SM field strength tensors. 

£ = ^ Xl\B^ + ||x7 5 xW^Wr + 0X7 *xG% v G<? (10) 

These operators are not velocity suppressed and they may naturally arise if the dark matter couple to the SM through 
loops of heavy messengers which carry SM charge HD. Unlike the previous model discussed, gauge invariance ensures 
that two of the above operators contribute to multiple annihilation channels. These two operators also have interfering 
contributions to the same annihilation channels. The coefficients fti/A 3 and /A 3 control the coupling of dark matter 
to four pairs of electroweak gauge boson final states: ZZ, Z 7 , 77 , and W + W~, see the relevant diagrams in Fig. [ 8 ] 



FIG. 8: Dark matter annihilation into gauge bosons. 

The thermally averaged annihilation cross section to di-boson final state is given in Eq. [TT] below. 


(ii) 


We now explore the bounds on the parameter space of this model, the space consisting of the dark matter mass m x 
and the three operator coefficients. For any given point in parameter space, the total photon flux due to annihilations 
results from the sum the annihilations into each of the five di-boson channels. We see that the operator with coefficient 
k 3 /A 3 controls the annihilation to gluons and only gluons, and therefore factorizes in the total annihilation rate as 
per the discussion before. However the total rate to other di-boson pairs depends on two operator coefficients. In this 
case, the total annihilation rate has the form (fc* = Kij A 3 ) 
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(akl + bkl) 2 + ck 2 = (crv) tot ( 12 ) 

We may again study the system by fixing the total annihilation rate as in the previous section. Unlike the case 
in the previous section, the rate constraint is not linear; however for each specific dark matter mass, the couplings 
which saturate the annihilation rate are still constrained to he on a 2-dimensional surface. This surface is a triangular 
section of an ellipsoid. The radii of the ellipsoid vary as the total annihilation rate is changed. We can again present 
dwarf constrains on the space of operator coefficients using a simple visualization, now on an ellipsoidal surface rather 
than a triangle. In Fig 9, we show the ellipsoidal surfaces which satisfy the constraint that the total DM annihilation 
rate is 10 times the thermal rate. 

We consider four different dark matter masses for which this constraint is satisfied, m x = 10, 75,150 and, 200 GeV. 
For each DM mass, certain combinations of the operator coefficients will satisfy the rate constraint. As the DM mass 
is decreased, lower cut-off scales A will obey the annihilation rate constraint, just as in the previous section. For each 
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point on the ellipsoid, we may check if the combined photon flux from di-boson annihilations violates or satisfies the 
dwarf constraints. Just as before, we find that parameter regions of low mass (hence low effective cut-off scale) are 
ruled out. In our scan of Dark Matter masses, we found that there is no allowed parameter space for DM masses 
much below 10 GeV which also satisfies the constraint (crv} tot = 10 {cfv)th- 





FIG. 9: Limits on effective operator coefficients for 10 times the thermal annihilation rate given by Fermi dwarf bounds. The 
axes show the value of fc/A 3 . Bounds on operators are shown for four dark matter masses, m x = 10,75,150 and, 200 GeV 
(top left, top right, bottom left and bottom right respectively). Red points indicate that the operator coefficient is ruled out, 
blue points are unconstrained. 

In the Fig. [9] plots, the red points are ruled out by the Fermi-LAT data, while the blue points are unconstrained. 
As we scan through parameter space increasing the Dark Matter mass different annihilation channels open up which 
contribute to the total photon flux. The first plot gives constraints for a 10 GeV dark matter particle. At such low 
masses m x < mz/2, only annihilations to gluons and photons are kinematically accessible. The photon energy E 1 is 
equal to the DM mass m x and the operator coefficient values that yield a high di-photon annihilation rate and are 
ruled out. For DM masses above half of the Z mass, the Zy annihilation channel open up and starts to increase in 
magnitude but we see the the characteristic photon energy of the Zy annihilation is less than that of the photons 

in the yy annihilation channel, here E 1 = m x ^J 1 — 4m^/ m 2 z . Therefore the electroweak operator coefficients in 
this mass region are not too strongly constrained. For higher dark matter masses, m x > mw/z-, the WW and ZZ 
annihilation channels also open up as we see in the bottom two plots where fi x is 150 and 200 GeV respectively. 

We have so far proceeded by fixing the total annihilation rate. We will now present the same constraints in a semi- 
orthogonal plane of parameter space where the total DM annihilation rate is allowed to vary. In Figs. [l 0 ]and[TT] we 
show the Fermi-LAT bounds in plane of operator coefficients. For simplicity in the figures we made the redefinition, 
ki = Kij A 3 . We now allow the coefficients k\ = fiq/A 3 and k 2 = «$/A 3 to vary freely. Here, ks is determined at 
each point by fixing the fractional annihilation rate to gluons R gg . We may then scan through increasing masses and 
determine which points in parameter space are ruled out. Note, we are not holding k% = k 3 /A 3 constant; as m x is 
varied this parameter must be varied in order for R gg to remain at a fixed value. Each point in the plot corresponds 
to a unique value of the three operator coefficients so each specifies a particular total annihilation rate for DM. As in 
the previous section, the parameter space exclusions are derived from the upper bounds on y ray flux from dwarfs. 
For fixed operator coefficient values, the annihilation cross section goes as m^, so the total photon flux for each point 
increases goes like m^. We may thus put an upper bound on DM masses for fixed values of operator coefficient. 
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We plot upper limits on Dark Matter masses in Fig. [lO] The top left, top right, and bottom plots in Figs. [10| 
and 11 correspond to 0%, 30%, and 70% annihilation rate into di-gluons respectively. The solid, dashed, and dotted 
lines represent the regions where Rzz, and i7 77 go to zero respectively. The vanishing of each of these channels 
coincide with a specific choice of the couplings in Eq. [To| given by 


^ 1 =^ 2 =^ Rz'y = 0 
s 2 

k 2 = —Rzz = 0 

c l 

c 2 

hv2 = =>• i ^ 77 = 0 (13) 

S w 

As one would expect, the 77 channel tends to be the most powerful in constraining annihilation rates when allowed. 
As we can see in Fig. [lOj much of the parameter space will be ruled out by the 77 channel since the limits are 
stringent enough that other channels do not turn on before photon lines rule them out. 



FIG. 10: These c ont ours correspond to upper bounds on mn x (in GeV) from annihilation into all final states allowed from the 
operators in Eq. 10 with ki = mf A 3 . The amount to the gg final state is fixed in each of these plots as a percentage of the 


total rate, 0% (upper left), 30% (upper right), 70% (lower). The solid, dashed, and dotted lines correspond to a combination 
of ki and K 2 which turn off the final states Z 7 , ZZ , and 77 respectively. 
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FIG. 11: These contours correspond to upper bounds on (crv)tot in units of the thermal relic, (av) Th, from annihilation into 
all final states allowed from the operators in Eq. [lojwith ki = m/ A 3 . The amount to the gg final state is fixed in each of these 
plots as a percentage of the total rate, 0% (upper left), 30% (upper right), 70% (lower). The solid, dashed, and dotted lines 
correspond to a combination of k\ and k 2 which turn off the final states Z 7 , ZZ , and 77 respectively. 


Notice that when the operator coefficients are very small, flux bounds allow for much larger DM masses. If we look 
at the region where the 77 final state is off in parameter space, we still expect annihilations into the final states with 
WW/ZZ and Z 7 which allow for strong constraints to be placed. In Fig. 

DM annihilation rate in the plane of effective operator coefficients. 
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we plot the upper bound on the total 


Limits of the EFT and Collider Constraints 

We will now say a few words about the maximum range of validity of the EFT for DM annihilation processes. The 
most natural completion of a model with a vector boson portal involves loops of heavy messenger particles which are 
charged under the SM gauge groups such as those found such as in m- For the EFT to be valid, the momentum 
transfer of the process should be less than twice the messenger mass. Since the square root of the c.o.m energy of the 
process is twice the Dark Matter mass, we get m x < M mess . We assume that the loops roughly scale like f&ig x /M^ ess 
where f is some numerical factor, M mess a messenger mass, and g x the coupling of the messenger to the dark matter. 
Therefore, our effective operator coefficient ki/A 3 corresponds to /c^g 2 /M^ ess . We can then estimate the messenger 
masses in terms of our effective cut-off for the maximum perturbative value of the hidden sector coupling g x . We find 

Mmess A (/aM 67 r 2 ) . For EFT validity m x < A (/aM 67 r 2 ) 1 ^ 3 . For order 1 values of f, the validity limits for 

the U(l) operator for example are roughly m x < A; in regions of parameter space where the DM mass larger than 
the effective cut-off, derived bounds are not reliable. 

We now wish to synthesize our parameter space constraints with EFT-validity limits and with constraints from 
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colliders. We will first reinterpret our results given in Fig. [IT] in the plane of dark matter mass vs operator suppression 
scale for several specific combinations of operator coefficients and ^ 3 . First, we will assume ^3 = 0, that is there 

are no annihilation of DM into gluons. We now have a three dimensional parameter space of ra x , k\ /A 3 , /A 3 . Our 
exclusions in this parameter space from Fermi-LAT data are shown in Fig. [13] for certain simple ratios of couplings 
Ki and Kj 2 • The excluded regions of parameter space lay under the colored lines. The dot and dashed lines in Fig. [13] 
correspond to curves of constant total annihilation rate. 
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FIG. 12: Constraints from Fig. |11| reconfigured for the effective cut-off vs dark matter mass plane. The regions below the solid 
lines are excluded. The dashed lines correspond to curves of constant (crv)to t 


We recall that for the operators in Eq. 10 ( av ) oc rad /A 6 and thus upper bounds on annihilation cross section or 


dark matter mass correspond to lower bounds on the suppression scale A. Below the solid purple line is the region of 
breakdown of EFT validity for the indirect annihilation process for the U(l) operator. Below the solid brown line is 
the region of breakdown of EFT validity for the indirect annihilation process for the SU(2) operator. 

Vector boson portal models like those in this section have been widely studied in the context of collider production. 
At LHC, operators of this type will lead to the production of dark matter particle pairs in association with a single 
vector boson, pp XX + V where V is a gluon, photon, W or Z. No mono-boson signal is currently in excess and 
therefore collider constraints place lower bounds on the effective operator cut-offs A * for each given dark matter mass. 
Studies on mono-boson collider signatures for the pseudo-scalar DM current may be found, for example, in references 
[ 22 ] and [ 2 T]. In the plot below we consider the effects of the SU( 2 ) operator only for a simple comparison of collider 
and indirect detection bounds. We show our exclusion contour from combined dwarf data, along with the ATLAS 
constrain from the mono-W channel (given by [22]) for the pseudo-scalar current operator. We have also plotted 
contours of constant total DM annihilation rate. 


Limits on A with k 1 =1, k 2 = 0 

A(GeV) 



— Collider mono-W channel 

— Combined dwarf limit 

— Max EFT validity k^= 0, x 2 =max 

0.1 <C7V> 

. <OV> 

- 5<ov> 

— — 10<ov> 

- 100<ov> 


FIG. 13: Constraints from Fig. [IT] reconfigured for the effective cut-off vs Dark Matter mass plane 
We see that following lines of constant annihilation rate we arrive at a minimum dark matter mass bound. Following 
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lines of constant operator coefficient we arrive at a maximum DM mass bound. The Fermi-LAT dwarf bounds are 
more constraining than collider analyses for DM masses above a few hundred GeV. It may be true that these bounds 
are more strict than collider bounds for regions of small A as well. The validity of the EFT prescription for the collider 
bound requires that the momentum transfer yfs < 2M mess . The general procedure for collider bounds is to truncate 
the number of events, eliminating those where yjs is larger than 2 times the maximum messenger mass [8]. Previous 
collider analyses for models with Dll type gluon operators estimate that the region of maximum EFT validity limits 
interpretation for cut-off scales below 350 GeV. For these low effective cut-off scales, event rates are severely truncated 
[23] . In our analysis for the SU(2) and U(l) operators, cut-off scales of around 500 GeV, imply messenger masses 
at or under roughly 1 TeV. We generated a sample of mono-W events for the process pp W\X a t 8 TeV c.o.m. 
energy using in MADGRAPH [24] . Barring any other cuts, we found that less than 10 percent of the events pass the 
above truncation requirement for an effective cut-off of 500 GeV and DM masses of 100 GeV. We do not know how 
far under 500 GeV the collider limits hold for these operators and this is a matter for further investigation. 


CONCLUSIONS 

We have given indirect detection constraints derived from Fermi-LATs dwarf spheroidal galaxy data on variety of 
models with multiple final state annihilations of DM. These models included dark matter portals from third generation 
fermions, SM vector bosons, and popular effective field theory models formally analyzed in collider physics. We have 
presented constraints in model space on surfaces of fixed total annihilation rate, and in planes of freely varying effective 
operator coefficients. In planes of fixed annihilation rate our bounds take the form of lower mass limits for dark matter 
mass/effective operator cut-off. Our mass bounds are in the 10-100s of GeV range, and we have presented these results 
using a visualization on triangles and ellipsoidal sections. We have also compared our constraints to those obtained 
by collider production for vector boson portal models. In general we find that dwarf constraints overtake collider 
constraints for dark matter masses greater than a few hundred GeV, and in regions of low effective cut-off where 
there is theory break down for collider applications. While existing results indicate that EFT validity fails under a 
few hundred GeV in collider studies, the exact range of maximum validity in vector boson portal models is yet to be 
determined. This is a course for future work. 

Dwarf constraints used are among the most limiting constraints for Dark Matter models, and therefore present 
the best possibility of ruling out regions of Dark Matter parameter space. We reiterate that a more sophisticated 
combination of dwarf data, which doesn’t assume annihilation into only one channel would provide much tighter 
constraints on our parameter space, and this is another avenue for further work. We also give another caveat here, 
we have implicitly assumed the most uncontroversial dark matter profiles for dwarf galaxies. Assumptions of different 
Dark Matter profiles will lead to wiggle room in our constraints and require further study. 


APPENDIX 


We present lower limits on the effective cut-offs A& and A r for our third generation fermionic model. We present 
lower limits for the fixed total DM annihilation rate 10 (crv)Th- 
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FIG. 14: Lower bound on operator coefficients and A* 
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